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ABSTRACT 

Recent discoveries of recombining plasmas (RPs) in supernova remnants (SNRs) have dramatically 
changed our understanding of SNR evolution. To date, a dozen of RP SNRs have been identified 
in the Galaxy. Here we present Suzaku deep observations of four SNRs in the Large Magellanic 
Cloud (LMC), N49, N49B, N23, and DEM L71, for accurate determination of their plasma state. 
Our uniform analysis reveals that only N49 is in the recombining state among them, which is the 
first robust discovery of a RP from an extra-galactic SNR. Given that RPs have been identified only 
in core-collapse SNRs, our result strongly suggests a massive star origin of this SNR. On the other 
hand, no clear evidence for a RP is confirmed in N23, from which detection of recombination lines 
and continua was previously claimed. Comparing the physical properties of the RP SNRs identified 
so far, we find that all of them are categorized into the “mixed-morphology” class and interacting 
with surrounding molecular clouds. This might be a key to solve formation mechanisms of the RPs. 
Subject headings: ISM: abundances - ISM: individual (DEM L71, N23, N49, N49B) — supernova 
remnants — X-rays: ISM 

1. INTRODUCTION 

X-ray observations of supernova remnants (SNRs) al¬ 
low us to make accurate measurements of plasma condi¬ 
tions and elemental abundances in the supernova (SN) 
ejecta, providing unique insights into the progenitor’s 
evolution and explosion as well as the dynamical evo¬ 
lution of the remnants themselves. The hot plasmas 
in X-ray-emitting SNRs are commonly in a state of 
non-equilibrium ionization (NEI), where the ionization 
degrees of heavy elements are inconsistent with those 
expected for a collisional ionization equilibrium (CIE) 
plasm a with a certain electron temperature (e.g., iMasail 
Il984l ). It had been widely believed for a long time 
that the SNR plasma is always in low-ionization, and 
slowly ionizing to achieve the CIE (hereafter “ionizing 
plasma”; IP). In fact, recent X-ray observatio ns with sen¬ 
sitive satellites, like Suzaku (iMitsuda et ali 12007 1. have 
confirmed that the immediate post-shock gas in young 
SNRs indeed consists of extremely low ion ized atoms 
together with hot elec trons (e.g., lUchida et all 120131 
lYamaeuchi et al.ll2014bll . 

However, earlier ASCA observations had suggested 
presence of a recombining plasma (hereafter RP, where 
the atoms are overionized compared to the observed elec¬ 
tron tem perature) in a couple of SNRs, IC 443 and 
W49B dKawasaki et al.1 120021 . 120051) . The conclusive ev¬ 
idence for the RPs was revealed by later Suzaku obser¬ 
vations that discovered enhanced radiative recombina- 
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tion c ontinua (RRCs) in the X-ray spectra of the se rem- 
nants (j Yamaguchi et al .1120091 [Ozawa et al .1120091 ). These 
results were followed by observational studies of other 
SNRs, leading to a significant increase in the number of 
RP SNRs. To date, RPs have been discovered from a 
dozen of SNRs, which means that the presence of RPs is 
no longer unusual or surprising. 

The well-identified RP SNRs are all categorized into 
the so-called mixed-morphology (MM) class, defined 
by centrally-filled therma l X-ray e mission with a syn¬ 
chrotron radio shell (|Rho fe Petrel IT998h . While more 
than 25% of the X-ray- detected Gala ctic SNRs are clas¬ 
sified into this type (l.Tones et al.lll998f l , a physical process 
to form such characteristic morphology is still unclear. It 
is theoretically pointed out that dense ambient materi¬ 
als play an impo rtant role in the dynamical evoluti ons of 
MM SNRs fe.g.. lWhite fe Lorndl 199 ll lPetrukll200 111 . No¬ 
tably, there is another prediction that the formation of 
RP can also be explained by interaction between the SN 
ejecta and d ense mate rials surrounding the progenitor 
(e.g.. lItoh fc MasaM989 ). Therefore, physical processes 
that create the mixed-morphology and RP seem related 
with each other. 

It is noteworthy that the most RP SNRs identified so 
far were previously consi dered to have an IP or a nearly- 
CIE plasma (e.g., W44: ITTchida et all 12012 1. This im¬ 
plies that there are still a number of SNRs of which 
plasma state has been misclassified. Here, we present 
uniform analysis of high-quality data of MM SNRs (as 
well as typical shell-like SNRs simultaneously observed) 
in the Large Magellanic Cloud (L MC) obtained by the 
X-ray Imaging Spectrometer (XIS: fKovama et aD i 2007 ) 
on board Suzaku , in order to search for a RP in these 
SNRs. The LMC is particularly suitable for such sys- 
tematic stud y because of its low foreground extinction 
(iDickev fe Lock manl I1990D and its known distance of 
5(1 knc TPhastlll 99911 : we use this distance throughout the 
paper. 
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In §2, we outline our observations and data reduction 
procedures. The results of spectral analysis of four LMC 
SNRs, N49, N23, N49B, and DEM L71 are presented 
in §3 — where we will show a new discovery of a RP 
from N49. We discuss our results in §4, and finally give 
conclusions and future prospect in §5. The errors quoted 
in text, tables, and figures are at the 90% confidence level 
unless otherwise stated. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed two regions in the LMC aiming at the MM 
SNRs N23 and N49 during the Suzaku Cycle 8 phase. 
Detailed information of the observations is summarized 
in Table Q] In addition to the targeted sources, nearby 
SNRs DEM L71 and N49B are detected in each field of 
view (FoV) of the XIS. We analyzed their spectra as well, 
for comparison of the plasma states among the targets. 
Both DEM L71 and N49B are c l assified as standard shell¬ 
like SNRs dHuehes et al.l 120031 iPark et all l2003bf) from 
which a recombining plasma has never been observed to 
date. 

We analyzed data from three active CCDs; one is back- 
illuminated (XIS1) and the others are front-illuminated 
(XISO and XIS3), although only the merged and averaged 
front-illuminated (FI) CCD spectra are shown in Figures 
in the subsequent sections. We used HEAsoft tools ver¬ 
sion 6.11 for the data reduction. The uncleaned data 
were reprocessed using the calibration database released 
in September 2011, and screened with the standard event 
selection criteria for cleaned events. The resulting effec¬ 
tive exposures were 102 ks and 185 ks for the observations 
of N23 (with DEM L71) and N49 (with N49B), respec¬ 
tively. 

3. ANALYSIS AND RESULTS 

We show in Figure [T] the 0.6-10.0 keV images of the 
two observations after subtraction of non X-ray back¬ 
groun d (NXB) generated using xisnxbgen dTawa et al.1 
120081) . Since the diamet ers of the SNRs (< 1'.5; e.g., 
N49B: [Park et a0 12003b .) are smaller than the angular 
resolution of the X-Ray Telescopes of (XRTs; the half¬ 
power diameter ranges from 1'.8 to 2'.3), our analysis 
focuses on spectra from the entire SNRs. We extract the 
source spectra from the regions enclosed by the circles 
shown in Figure [l] The background spectra are taken 
from the entire FoV excluding the source regions as well 
as the CCD corners illuminated by the 55 Fe calibration 
sources. 

Figure [2] compares background-subtracted spectra of 
the four SNRs, showing that N49 is the brightest among 
them (with the count rate about ^10 times higher than 
the others). Owing to the good energy resolution of the 
XIS, we clearly detect K-shell emission lines of Ar, Ca, 
and Fe from this remnant, for the first time. The spec¬ 
trum of N49 exhibits strong Lyct lines of Mg, Si and S, 
whereas there is little or no signal of these lines from the 
other remnants. This indicates that the heavy elements 
in N49 are more highly ionized than those in the others. 
As our main goal is to identify a RP SNR(s), we first 
perform detailed spectral analysis of N49 (§3.1), which is 
followed by analysis of the other remnants (§3.2-3.4). 

We use the S PEX software version 2.04.01 
dKaastra et al.l 179961) for the spectral fitting by taking 



Figure 1. Vignetting-corrected XIS FI images of the N23 (left) 
and N49 (right) regions in the 0.6-10.0 keV band with an equatorial 
coordinate grid. DEM L71 and N49B are also detected in the FoV 
of the XISs indicated by the solid squares. 



Energy (keV) 

Figure 2. XIS FI spectra of N49, N23, N49B and DEM L71, 
where the background spectra are subtracted. The spectra of N49B 
and DEM L71 are multiplied by factors of 0.1 and 0.01, respec¬ 
tively. The arrows represent the line centroids of Mg-Lya, Si-Lya 
and S-Lya. 

into account the detector and telescope responses (so- 
called “redistribution matrix” and “ancillary response”, 
respective ly) g enerat ed by xisrmfgen and xissimarfgen 
(jlshisaki et al.l 120071) . The data around the neutral 
Si K-shell edge (1.77-1.83keV) are ignored because of 
the poor accuracy in the response function at these 
energies 5 . During the analysis, we separately set ab¬ 
sorption column densities in the Milky Way (-/V H (mw)) 
and the LMC (-/V H ( lmc))- The former value is fixed to 
6 x 10 20 cm -2 (iDickev fe Lockmanl 1199(31) . Elemental 
abundances of the latter component a re fi xe d to the av¬ 
erage s values of the LMC (~0.3 solar; IRussell fe Dopital 
fl992) . 

3.1. N49 

N49 is one of the brightest SNRs in the LMC at the var¬ 
ious wavelengths (e.g.. lLong et aljfl981l : IDickel fe Milnel 
1998). The radio con tinuum and Hct emission exhibit 
a clea r shell structure dVancura et al.l179921 : IDickel et al.l 
1995), whereas the X-ray image shows an irregular mor¬ 
phology brightest inside the radio shell in the southwest 

5 http://heasarc.nasa.gov/docs/Suzaku/analysis/sical.html 
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Table 1 

Observation logs. 


Target 

Obs. ID 

Obs. Date 

(R.A., Decl.)j 2 ooo 

Exposure 

N23 & DEML71 

807008010 

2012 April 4 

(76.45, -67.96) 

102 ks 

N49 & N49B 

807007010 

2012 May 9 

(81.50, -66.08) 

185 ks 


N49 
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Figure 3. XIS FI spectrum of N49 (top panel). The best-fit 
model is overlaid with the black solid line. The solid blue and red 
lines represent the best-fit high-/cT e (RP) and the low -kT e (IP) 
components, respectively. The orange and green lines show the 
power-law and the blackbody components for SGR 0526—66, re¬ 
spectively. The dotted lines show Gaussians for the Fe L emission 
missing from the plasma code we used. The residuals are shown in 
panel-(d). Panels (a), (b), and (c) represent the residuals from the 
models of single IP, two IP, and two IP plus Gaussians, respectively. 


region dPark et al.li2012f) . which is the typical character¬ 
istic of a MM SNR. N49 is spatially overlapping with the 
soft gamma-ray repeater (SGR) 0526—66 (jCline et all 
1982), although the physical ass ociation between the 
SGR and SNR is u nder debate (jGaensler et al.l 120011: 
iBadenes et al.ll2009l) . Using Chandra data. iPark et al.1 
(|2012f) revealed that the X-ray spectrum of SGR 0526—66 
is best reproduced by a model consisting of a blackbody 
(BB) with kT BB = 0.44 ±0.02 keV and a power law (PL) 
with T = 2.50^0^2- Unfortunately, our Suzaku data do 
not spatially resolve the SGR from t he S NR. We thus 
add the best-fit model bv IPark et al.l ( 2012 1 to the model 
spectrum of N49, but allow its spectral parameters to 
vary within the reported statistical uncertainty. Since 
IPark et all (120121) found a significant long-term variabil¬ 
ity in the X-ray flux of the SGR, we also allow the fluxes 
of both BB and PL components to vary freely. 

We first apply a single-component IP model where the 
initial ionization balance is dominated by the neutral 
state for all heavy elements. The free parameters are 
the electron temperature ( kT e ), ionization timescale ( n e t , 
where n e and t are the electron number density and the 
elapsed time since the gas was shock heated), emission 
measure (EM = f n e n B V), and abundances of O, Ne, 
Mg, Si , S, Ar, Ca and F e relative to the solar values of 
lAnders fe Grevessel (U989D . The Ni abundance is linked 
to the value of Fe, while the abundances o f the other ele- 
ments are fixed to the LMC averages of lRussell fe Dopital 


(119921) . This model yields the best-fit electron tempera¬ 
ture of ~ 0.6 keV, but fails to reproduce the overall spec¬ 
trum with an unacceptable y 2 /d.o.f. of 4936/820. The 
residual is shown in panel (a) of Figure [3] The particu¬ 
larly large disagreement between the data and model is 
seen below ~1 keV, implying presence of another plasma 
component with a lower temperature. Assuming that 
this component originates from a swept-up interstellar 
medium (ISM) t hat is dominant at the o utermost region 
of the remnant (jPark et, al.l l2003al 120121 ) , we added an¬ 
other IP component with the abundances fixed to the 
LMC mean values. Although this two-component model 
gives a slightly better fit (y 2 /d.o.f. = 4297/817), some 
line-like residuals still remain as shown in panel (b) of 
Figure [3] The residual at ~0.8 keV, also observed i n our 
previous studies of othe r objects (e.g.. lllchida et all2013t 
iNakashima et all 120131 ). is likely due to the well-known 
uncertainty in the emissivity ratio between the Fe XVIII 
L-she ll emission of 3s —> 2 p and 3d — > 2 p (IGu et al.1 
120071 ). On the other hand, the residuals at ~ 1.2 keV 
and ~ 1.3 keV are caused by L-shell emission from high 
quantu m numbers (n > 5) mis sing from the plasma code 
SPEX dBrickhouse et al.ll2000D . We compensate for these 
lines by adding three Gaussians at these energies, and 
obtain the improved fit (x 2 /d.o.f. = 2366/815) with the 
best-fit kT e of ~ 0.6 keV and n e t of ~ 10 13 cm -3 s for the 
high -kT e component. The obtained ionization timescale 
indicates that this IP model actually represents a CIE 
plasma at the given temperature. 

The fit left the largest residuals at 2.0 keV and 2.6 keV 
(see panel (c) of Figure (3]), which correspond to the cen¬ 
troid energies of the Si Lyct and S Lyct emissions, respec¬ 
tively. This suggests that the average charge of these 
elements is higher than that expected in a 0.6-keV CIE 
plasma. Furthermore, a hump-like feature is found in 
the residual around 2.7 keV, indicating enhancement of 
the RRC of Si, as was observed in the other RP SNRs 
(e.g.. IC443: lYamaguchi et a.1.1 l2009h . Therefore, we in¬ 
troduce a RP model (using an “NEIJ” model in SPEX 6 ) 
for the high-fcTe component. In addition to the parame¬ 
ters given in the IP model, the RP model has another free 
parameter kT m { t of higher than kT e , which describes a re¬ 
combining history by n e t in the kT e plasma starting from 
the initial ionization temperature of kT lni t. This model 
dramatically improves the fit (\ 2 /d.o.f. = 1270/813), 
removing the large residuals above 2 keV as shown in 
the panel (d) of Figure 02 The best-fit model compo¬ 
nents and parameters are given in Figure [3] and Tabled 
respectively. 

Interestingly, the RP model reproduces the Fe Ka 
emission at ~ 6.6 keV as well, although free electrons 
in the 0.6-keV plasma are not energetic enough to excite 

6 See http://www.sron.nl/files/HEA/SPEX/manuals/manual.pdf 
for details. 
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Figure 4. XIS FI spectra of N23 (upper left), N49B (upper right) and DEM L71 (bottom). The best-fit models are overlaid with the 
black solid lines. The solid blue and red lines represent the best-fit high-fcT e and the low -kT e components, respectively. The dotted lines 
show Gaussians for the Fe L emission missing from the plasma code we used. The orange line in the model of N23 shows the power-law 
component for CXOU J050552.3—680141 (see text). 


K-shell electrons of Fe. This indicates that the observed 
Fe Kct emission predominantly originates from cascade 
processes after the radiative recombination into the ex¬ 
cited levels of He-likc ions (i.e., Fe 25+ + e~ —> Fe 24+ *). 
We emphasize that this is why kT- m it is constrained to the 
very high value (> 10 keV) so that a significant fraction 
of the Fe ions remains at the H-like state in the cur¬ 
rent plasma. It should also be noted that the observed 
centroid energy 6629+26 eV is consistent with the value 
expected from the best-fit RP model (6658 eV) but sig¬ 
nificantly lower than that for a typical CIE plasma at 
kT e = 5-10keV (~ 6680 eV). This is another piece of 
evidence for the RP; the forbidden and intercombination 
lines are likely to be enhanced by the recombination pro¬ 
cesses. 

3.2. N23 

In the pr evious work us i ng XM M-Newton MOS and 
RGS data, iBroersen et all (2011 ) claimed that the X- 
ray spectrum of N23 can be best reproduced by a two- 
component NEI model of which one of the components 
was a RP with kT m i t = 3.0 keV and kT e = 0.18keV. We 
first fit the XIS spectrum fixing their best-fit parame¬ 
ters, but fail to reproduce the XIS spectrum (y 2 /d.o.f. 
= 455/236), leaving large residuals around the Ne-Lya 
emission and the Fe XVII recombination edge at ~ 
1.25 keV. The fit is not improved even if we allow the 
abundances to vary. Therefore, we thaw the tempera¬ 
tures and ionization parameters (but restricting kT- m \t 
to be higher than fcT e ), obtaining a significantly better 
fit with y 2 /d.o.f. = 224/234 and n e t > 10 12 cm -3 s for 
both components. This high n e t value indicates that the 
plasmas are in nearly CIE independently from the initial 


ionization population (otherwise ionizing). In fact, we 
are not able to constrain kT- m it with this model. Finally, 
we apply an IP model consisting of a LMC-abundance 
component (for ISM) and a free-abundance component 
(for ejecta), and obtain a slightly better fit (\ 2 /d.o.f. 
= 218/234) than the RP (or CIE) model. The best-fit 
parameters and models are given in Table [2] and Fig¬ 
ured! respectively. During the analysis, we fix the spec¬ 
tral shape and flux of the faint X-ray compact source 
CXOU J050552.3—680141 associated with N23 to the re¬ 
ported values of lHavato et al.l (|2006l ). although its con¬ 
tribution is negligible at <3keV (see Figure |4j- 

3.3. N49B 

N49B is a typical shell-like SNR containing no X- 
ray po int-like source within the SNR shell (Park et all 
(2003 5). We fit the XIS spectrum with the same two- 
component IP model applied to N49, obtaining the best- 
fit results given in Table [2j The n e t value for the ejecta 
component (~ 8 x 10 10 cm -3 s) is the lowest among the 
SNRs studied in this work. We reveal that the abundance 
of Mg is significantly higher than those of the other ele¬ 
ments (e.g ., Mg/Q) ~ 2.0, Mg/Ne 2.8, Mg/Si ~ 3.6). 

I Park et all (20031)1 ) previously found “Mg-rich” ejecta 
near the center of the remnant. Our result shows that 
the Mg enrichment is confirmed even in the integrated 
spectrum from the entire SNR. We find no evidence for 
a RP; if kTi n > kT e is assumed, the model clearly fail 
to reproduce the spectrum with x 2 /d.o.f. = 1582/426. 
Therefore, we conclude that this SNR is dominated by 
IPs. 


3.4. DEM L71 
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Table 2 

Best fit parameters 


Component 

Parameter 









N49 

N23 

N49B 

DEM L71 

Absorption 

-Vh(lmc) (xlO 21 

cm 2 ) 

q ci+0.01 

°- ol -0.04 

2.13 ±0.05 

2 30 +0 01 

^.ou_0 06 

q '7/~J + 0.30 

'-0.50 


-Vh(mw) (xlO 21 

cm -2 ) 

0.6 (fixed) 

0.6 (fixed) 

0.6 (fixed) 

0.6 (fixed) 

ISM 

kT e (keV) 


0.30 ±0.01 

0.17 ±0.01 

0.21 ±0.01 

0.16 ±0.01 


n e t (xlO 11 cm -3 

s) 

q 4q+0-09 

o.^y—0 26 

> 10 

-1 OqPO.OS 
±,zo —0.07 

> 10 


EM (xlO 58 cm - 

3 ) 

q 34+ 003 
o.O4_ 0 01 

17.6 ±0.3 

6.0 ±0.9 

I9/4 

Ejecta 

kT e (keV) 


0.62 ±0.01 

0.61 ±0.01 

0.79 ±0.01 

0.69 ±0.01 


fcT init (keV) 


11 ± 1 

0.01 (fixed) 

0.01 (fixed) 

0.01 (fixed) 


n e t (xlO 11 cm -3 

s) 

7 00 +0 - 44 

i -UU_Q 02 

5-08tg1? 

0.84 ±0.01 

1.77 ±0.07 


O 


1 ro+0.04 
x - oz, -0.03 

3.5 ±0.2 

0 63+u- 03 
U.DO — 0 02 

in ^ 

q q 
6 6 
+ 1 
O 
CO 

0 


Ne 


0-98t°;° 3 

1.25 ±0.07 

0.45 ±0.01 

0.25 ±0.02 


Mg 


0.82 ±0.02 

1.16 ±0.07 

1.24 ±0.02 

0.62 ±0.03 


Si 


-1 oq+0.05 

J-.oZ—o 02 

1.2 ±0.4 

0.34 ±0.02 

0.48 ±0.04 


S 


1 4f)+ 0,03 

1.25 ±0.08 

0.50 ±0.08 

1.0 ±0.2 


Ar 


1 o+O.l 
llCS -0.5 

< 3.1 

0.7 ±0.6 

1 c+3.1 
±,o -0.3 


Ca 


1 c+0.2 
1,o -0.9 

(=Ar) 

(=Ar) 

(=Ar) 


Fe (=Ni) 


0.24 ±0.01 

0.46 ±0.02 

0.35 ±0.01 

0.62 ±0.01 


EM (xlO 58 cm - 

3 ) 

1.73 ±0.01 

0.62 ±0.01 

0.79 ±0.01 

1.36 ±0.01 


X 2 / d.o.f. 


1270/813 

218/234 

659/423 

291/229 


DEM L71 is a middle-aged shell-like SNR, where 
a double-sho ck morp holog y was observed by Chandra 
observations (IHughes et al.l 120031) . suggesting the pres¬ 
ence of both revers e-shocked ejecta and swept-up ISM. 
iHughes et al.1 ( 2003 :1 also revealed that the X-ray spectra 
from the interior regions are dominated by strong Fe L 
emission. Similarly to the spectral analysis of N49B, we 
fit the XIS spectrum of DEM L71 with a two-component 
IP model consisting of a low-fcT e ISM and a high -kT e 
ejecta. The fit is acceptable with % 2 /d.o.f. of 291/229 as 
is given in Table [2j The Fe abundance is highest among 
the four SNRs, consisten t with the previous studies (e.g., 
IHughes et al.lll998Ll2003f l where the Type la origin of this 
SNR was suggested. If we assume fcTjmt > kT e , an un¬ 
acceptable fit (x 2 /d.o.f. = 399/228) is obtained, ruling 
out a RP scenario for this remnant. 

4. DISCUSSION 

4.1. Plasma Conditions in the Observed SNRs 

We have systematically analyzed X-ray spectra from 
the four LMC SNRs (N49, N23, N49B, and DEM L71), 
and revealed the robust evidence for overionization in 
N49 — this is the first discovery of a RP from an extra- 
galactic SNR. The other SNRs, including N23 from which 
prese nce of a RP was previously reported dBroersen et al.1 
Ml) , can be fairly well characterized by an IP or a 
nearly-CIE plasma. 

The previous claim of the recombining state in N23 
was based on the enhanced G-ratio of the O VII 
lines measured using the XMM-Newton RGS. Although 
iBroersen et al.l (l2011f) claimed that the broadband spec¬ 
trum of the XMM-Newton MOS can also be reproduced 
by a RP model, their best-fit model is clearly ruled out 
by our analysis of the XIS spectrum with better photon 
statistics. Since our analysis excluded the O VII emis¬ 
sion, it is still possible that the swept-up ISM (low-fcT e 
component) is partially recombining. We note, however, 
that there are several other processes that can enhance 
the G-ratio: resonance scattering and charge exchange. 


Notably, a similar G-ratio enhancement in the O VII 
lines was observed in DEM L71, and was indeed inter¬ 
preted to be a consequence of the resonance scattering 
(Ivan der Hevden et al.ll2003 ). Since the charge exchange 
proces s seems to work in s ome evolved SNRs (e.g., Pup- 
pis A: iKatsnda et al.ll2012fl . this process may also be re¬ 
sponsible for the high G-ratio observed in N23. 

4.2. Origin of the Recombining Plasma 

We summarize the characteristics of the four SNRs 
studied in this work in Tabic [3] Only N49 is i nteract¬ 
ing w ith dense cl ouds identified by t he CO dBanas et al.1 
119971) and Ha dMelnik fe CopettH 120131) observations. 
Given that similar cloud interaction is observed in most 
SNRs from which th e presence of a RP has been con¬ 
firmed (e.g., W49B; lOzawa et al.l 12009 ). thermal con¬ 
duction into surrounding clouds might play an impor¬ 
tant role in forming the RP (see I Zhou etahl l201ll for 
a discussion). In this scenario, there should be a spa¬ 
tial correlation between the electron temperature of the 
SNR plasma and the cloud density. Future deep observa¬ 
tions of N49 with better spatial resolution are necessary 
to assess this possibility. 

An alternative is that the adiabatic expansion of the 
SNR had caused a rapid cooling of electrons and the re¬ 
sulting recombining state of the plasma. In this scenario, 
the SN ejecta should interact with a dense circumstellar 
matter (CSM) in the early phase of the SNR evolution, 
and the highly ionized ejecta and CSM expand drasti¬ 
cally after the SNR blast wave breaks out to the low- 
density ISM (lltoh fe Masailll989l) . Therefore, the ioniza¬ 
tion timescale ( n e t) in this scenario is characterized by 
the elapsed time and the density evolution history since 
the break out took place. To evaluate the possibility 
of this scenario, we estimate the electron density in the 
present RP of N49 using the derived EM (= n e nnV). 
Since the plasma density is know n to be highly inhomo¬ 
geneous in this remnant dPark et al.1120121) . we analyze 
archival Chandra data of N49 to investigate the surface 
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Table 3 

Summary of the Four LMC SNRs Studied in This Work 


SNR 

RP 

Age (10 4 yr) 

SN Type 

Morphology 

Cloud Interaction 

References 

N49 

Yes 

0.4-0.5 

Core-collapse 

MM 

Yes (CO and Ha) 

this work, 1, 2, 3 

N23 

No 

0.4 

Core-collapse 

MM 

No 

2, 4 

N49B 

No 

1.09 

Core-collapse 

Shell 

No 

5 

DEM L71 

No 

0.4-0.5 

la 

Shell 

No 

6 , 7, 8 


References. — f I ') [Pa rk et al.||2012|; (2) Banag_et_alj[T997; (3) Melnik & Copctti 2013; (4) Broe rsen et ahlfeOllI ; (51 IPark et al.l l2003b; 
(6) IGhavamian et akll2003l ; (7) IHughes et ah 119981 ; (81 IHnghes et aljl2003[ 


brightness profile. We find that ~ 50% of the photon flux 
in the Si band (1.7-2.1 keV) is coining from the brightest 
southeast region (approximately a 0/4 x 0.'25 x 0.'25 el¬ 
lipsoid) , and the other ~ 50% is from the remaining faint 
region in the entire SNR (a sphere with a radius of 0.'6). 
Therefore, the volumes and the electron densities are es¬ 
timated to be 9.6 x 10 57 cm 3 and 1.0 cm -3 for the bright 
region, and 7.3 x 10 58 cm 3 and 0.38 cm -3 for the faint 
region, at a distance of 50kpc. If we simply divide the 
best-fit n e t in the ejecta component (7.0 x 10 11 cm -3 s) 
by the derived electron density, we obtain plasma ages 
of 22-58 kyr, more than a few ti mes l arger than the dy¬ 
namical age of N49 (~4800 yr; IPark et al]l2012i) . This 
is not surprising because the plasma must have h ad a 
higher electron density in the past. lYamaguchi et al.1 
(2012) also obtained a similar result from the Suzaku 
observation of IC 443, where the plasma age estimate 
(~ 11 kyr) is significan tly higher than the SNR age of 
~ 4000 yr (jTroia et al.ll2008D . 

In Table U we summarize physical properties of the 
RP SNRs identified so far. While the estimated ages are 
roughly correlated with fcT e , the values of n e t show no 
clear trend. This may be partly due to a technical reason 
that the parameter n e t is coupled to that of fcTinit in the 
spectral fitting. Nevertheless, Table [4] shows that the 
n e t values exceed ~ 10 11 cm -3 s in all SNRs, implying 
that the overionization did not start from the last few 
hundred years, but from the early epoch in the SNR. 
This fact likely favors the adiabatic cooling scenario for 
the most RP SNRs. 

4.3. Progenitor of Nf9 and SGR 0526-66 

The SN type of N49 has been controversial. IPark et ahl 
ll20T2h suggested a Type la origin based on the Si/S 
abundance ratio in the ejecta measured by the Chan¬ 
dra data. On the other hand, the environment of N49 , 
i.e., a nearby OB associat ion dChu feKennicuttl 119881) 
and young stellar clusters dKlose et al.l I2004T) . are more 
common in core-collapse SNRs. A recent systematic 
study of Fe K emissions in young and midd le-aged SNRs 
also ag rees on its core-collapse origin (|Yamaguchi et al.1 
l2014al) . Given that a RP has never been observed in a 
Type la SNR, our results strongly support a core-collapse 
scenario for this SNR. It should also be noted that the 
formation of a RP needs dense CSM, which is only made 
by a massive progenitor with a sig nificant mass loss rate 
(jMorival l2012t iShimizu et al.l I2012H . The total mass of 
the RP component is roughly estimated to be ~26M 0 
using the density and emitting volume derived in §4.2, 
which is in marginal agreement with the progenitor mass 
of 12.5-21.5 Mq constrained from the local star forma¬ 


tion history around N49 dBadenes et al.l[20091) . 

Since the core-collapse origin is strongly suggested for 
N49, the compact object SGR 0526—66, which is known 
to be a magnetar candidate (e.g., IPark et all 12012 1. 
might be associated with the SNR. The large distance 
between SGR 0526—66 and the geometric center of N49 
requires a high kick velocity (~ 1100 km s -1 ) for the 
SGR at the SNR age of ~ 4800 yr, and hence IPark et all 
( 2012 1 argued against the physical association between 
these objects. However, the theoretical prediction is that 
a magnetorotational core-collapse SN explodes asymmet¬ 
rically with bipolar jets and produces a ma gnetar wi th 
a large kick velocity up to ~ 1000 km s -1 (iSawai et al.1 
120081 ). Moreover, there are several compact objects lo¬ 
cated extremely off-center from the associated RP SNRs. 
For instance, IGR J11 014—6 103 is located outside the 
shell of MSH ll-61A. [Pavan et al.l (|2014f) indicated that 
this neutron star is escaping from MSH 11-61A at a ve¬ 
locity exceeding 1000 km s -1 . Also, PSR 1758—23 is 
known as a “runaway pul sar” which might be physi¬ 
cally associated with W28 ()Frail et, al.lll993fl . Two other 
compact objects CXOU J171419.8-383023 in CTB 37A 
(lAharonian et al.ll2008aff and XMMU J061804.3+222732 
in IC 443 (lOlbert et al.ll200lh . are also found far from the 
geometric center of their host SNRs. Taking these facts 
into consideration, we propose that SGR 0526—66 is the 
plausible candidate of the compact remnant of N49. 

5. CONCLUSIONS 

We have shown the clear evidence for the RP in the 
SNR N49, which is the first report of the robust de¬ 
tection of a RP from an extra-galactic SNR achieved 
by the high sensitivity and the good energy resolution 
of the XIS on board Suzaku. The other SNRs, includ¬ 
ing N23 where presence of a RP was previously claimed 
(Broersen et al. 2011), can be well characterized by an IP 
or a nearly-CIE plasma. A future observation of N49 us¬ 
ing high -resol utio n X- ray spectrometers, like ASTRO-H 
(iTakahashi et al]l2014D . is highly encouraged; our Suzaku 
result predicts the detection of the Fe RRC as well as the 
enhanced forbidden lines due to the recombination pro¬ 
cesses (see §3.1). We also emphasize the importance of a 
systematic study of the RP SNRs, which will help us un¬ 
derstand physical processes responsible for the formation 
of the RPs in more details. 

The authors thank to Dr. T. G. Tsuru for a careful 
reading of our manuscript. This work is supported by 
JSPS KAKENHI Grand Numbers 26800102 (H.U.) and 
2450229 (K.K.). 
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Table 4 

List of the RP SNRs in the Order of Increasing Electron Temperature 


SNR 

Age 

fc71 n it 

kT e 

n e t 

Cloud Interaction 

Compact Object 

References 


(10 4 yr) 

(keV) 

(keV) 

(10 11 cm -3 s) 




G359.1—0.5 

> 1 

0.87 + q'J“ 

0.29 ±0.02 

< 4.42 

Yes 


1. 2, 24 

G346.6-0.2 


5 (fixed) 

0 30”*"°' 03 
u.ou_o oi 

4 8”*" 0 ' 1 

Yes 


3, 25 

W28 

<4 

3 (fixed) 

o 4n+ 0 - 02 
u -^ u -0.03 

6.31 

Yes 

PSR 1758-23 

4, 5, 6, 24 

W44 

rsj 2 

1 07 +0 - 08 

1 -vi _o.06 

0.48 ±0.02 

6.76 ± 0.5 

Yes 

PSR B1853±01 

7, 8, 9, 24 

CTB 37A 

~ 1 

5 (fixed) 

0 4Q"*" 0 09 

u.^y_ 0 .o6 

13±? 

Yes 

CXOU J171419.8—383023 

10, 11, 12, 24 

3C391 


i -°-0.6 

0.495±0.015 

i4.ol 4 ;® 

Yes 


13, 24 

MSH 11-61A 

1-2 

5 (fixed) 

rv r-i o —|— 0.0U4 

U.O-LO_Q QQ3 

12.2 ±0.4 

Possible 

IGR J11014-6103 

14, 15, 16, 26 

N49 

0.4-0.5 

ii ± i 

0.62 ±0.01 

7 00 +0 ' 44 

1 .UU_Q Q2 

Yes 

SGR 0526-66 

this work, 17 

IC 443 

~ 0.4 

10 (fixed) 

0.65 ±0.04 

9.8 ±1.1 

Yes 

XMMU J061804.3±222732 

18, 19, 20, 21, 24 

W49B 

~ 0.4 


±.oz o 02 


Yes 


22, 23, 27 


Note. — The plasma parameters for IC 443 and W49B were determined only from a high-energy band around Fe-K lines (e.g., 
IQhnishi et al.|[2014l : IQzawa et al.H2OO90 . while those for the other SNRs were determined from intermediate-mass elements. 


References. — (1 ) lAha ronian et al. 2008 E|; (2) | Ohn ishi ct al. 20 1 1]; (3) lYa mauchi et al. 2013; (4) Rho & Borkowski 2002]; (5) 

Sawada & Koyama 2012]; (6) “r unaway” Puls ar; Frail ct al. 1993; (7) Cox ct al. 1999; (8) Uchida ct al. 2012]; (9) Pulsar plus Pulsar Wind 
Nebul a; I Petre e t all 12 0021; (T O) lYam auchi et al. 2014; (11) Wolszczan et al. 199ll ; (12) Pulsar Wind Nebula candidat e; lAharonian et all 
2008a; (13) Sato ct al. 2014; (14) Slanc ct al. 2002; (15) Kamitsukasa et al. 20 l]|; (16) Pulsar p lus Pulsar Wind Nebula: IPavan et al.1 12014: 
(17) p ark ct al. 2012; (18) Lee ct al. 2008; (19) Yamaguchi ct al. 2009; (20) Ohnishi ct al. 201-1; (21) Neutronstar plus Sync hrotron Nebula; 
Olbcrt ct al. 2001; (22) Hwang et al. 2 0001 ; f23) IQzawa et a l. 2009: (24) Hewitt ct al. 2008: (25) Korale skv et al .11 19981 ; (26) Filipovic et all 
20051: (27) Keohane et al .1120071. 
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